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Abstract 

We present a microscopic derivation of the equation of motion for a vortex 
in a superconductor. A coherent view on vortex dynamics is obtained, in 
which both hydrodynamics and the vortex core contribute to the forces acting 
on a vortex. The competition between these two provides an interpretation 
of the observed sign change in the Hall angle in superconductors with mean 
free path I of the order of the coherence length £ in terms of broken particle- 
hole symmetry, which is related to details of the microscopic mechanism of 
sup er conductivity. 
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In recent years the interest in vortex motion in superconductors has revived, mainly 
due to the advent of high temperature superconductors (HTSC). As a consequence of the 
peculiar material properties, the physics of vortices in HTSC shows many new aspects not 
encountered in conventional superconductors!. A major theme is the sign change in the 
Hall effect in the superconducting state, as is observed in both HTSCSI and conventional 
superconductors! for temperatures T just below Tq- This Hall-anomaly cannot be under- 
stood within the framework of the conventional Bardeen-Stephen@ or Nozieres-Vineni theo- 
ries for vortex motion, that predict the Hall effect in the superconducting and normal state 
to have the same sign for all temperatures. Several attempts at a theoretical understanding 
of the phenomenon have been undertaken!!!, but none of these seem to explain the experi- 
mental data. In Ref.l Hagen et ai, comparing a number of experiments, conclude that the 
sign change of the Hall effect is an intrinsic vortex property that occurs if the electron mean 
free path I is of the order of the coherence length £. Within a phenomenological analysis, 
Feigel'man et a/.! interpret the sign change in terms of broken particle-hole symmetry and 
obtain good agreement with the experimental signatures! of this effect. It is the purpose of 
this Letter to report on a microscopic calculation of the dynamical single vortex properties 
that yields a unifying description of the physics involved and puts the results of the analysis 
of Ref.l on a firm theoretical basis. 

Before presenting the microscopic theory, we discuss our main results for the vortex 
equation of motion and the resulting Hall force and angle. In general one expects the forces 
on a vortex to consist of two contributions, i. e., one from the electronic states in the vortex 
core and one from the hydrodynamic flow far away from the core. The vortex equation of 
motion has the form 

[M c + M H ]V + 7] C V = {k h V t - {-ic + 1h}V} x z . 

Here V denotes the velocity of the vortex and Vr is the transport velocity due to an applied 
current density j = %Vt/$ [we consider a film or layered structure, the extension to a 3D 
geometry is straightforward]. The equation of motion includes a vortex mass M, a damping 



term 77, and the Lorentz and Hall force coefficients k and 7. We made a clear separation 
into core and hydrodynamic contributions by writing subscripts C and H respectively. Ex- 
trinsic forces due to pinning and the interaction with other vortices add to the r.h.s. of the 
above equation of motion, however, here we consider only intrinsic vortex properties. The 
coefficient for the Lorentz force, k^j = nn s , arises from the hydrodynamic flow around the 
vortex with the superfluid density described by n s i. The mass terms were considered by 
Suhl in a Ginzburg-Landau approach0. The core contributions rjc and 7c were calculated 
by Kopnin et ct/.0, who found [we use % = c = ks = 1] 

where n e denotes the electronic density, u ~ A 2 /e^ is the level spacing between the localized 
Caroli-de Gennes-Matricon (CdGM) states in the cor and r r is the relaxation time. 

The key point in the determination of the Hall angle anaii (tananaii = l/v) is to find 
the hydrodynamic contribution 7^. For comparison we remind the similar procedure for 
uncharged Bosons like 4 He, where the vortex core has no internal structure and the Hall 
force arises from the first order time derivative ipidtip in the Lagrangian density^. With 
%p = y/ne lip ', n the mean particle density, the corresponding contribution to the Lagrangian 
is 5C = —mp. In the presence of a vortex at R, the phase configuration is ip(r, r) = 
(p v (r — R(t)) with ^(r) = arctan(?//x). The Euler-Lagrange equation yields a Hall force 
FHaii — — 27rnV x z, or 7^ = 2im for Bosons. If no normal fluid component is present at 
T = the Hall and Lorentz force combine into the Galilei invariant Magnus force Fm = 
k(V t - V) x zi. 

A hydrodynamic contribution to the Hall force in a superconductor arises also from a 
first order time derivative in the Lagrangian. This is most clearly seen in a time dependent 
Ginzburg Landau (TDGL) approach, where a term 5C = (N' e /2AN e )Aid t A appears in 
the Lagrangian density^. This term depends on the electronic band structure through the 
derivative of the density of states N e at the Fermi level N' e = d^N e {n) \ tt=ep and is thus related 
to particle-hole asymmetry. Here A is the strength of the attractive BCS model interaction. 



Note that in BCS theory 2N' e /(AN e ) = N e d^hiTc- The same procedure as for 4 He leads 
to a small hydrodynamic contribution 7^ = 7r(A^/AAf e )|A| 2 of order n e (A/e^) 2 (N' e ~ 
n e /e 2 F ). Its exact magnitude and sign depends on the (experimentally accessible) details of 
the electronic band structure. Although core physics is lacking in a TDGL approach^, we 
will see in the following that TDGL does predict the correct hydrodynamic contribution jjj. 
A hydrodynamic contribution to 7 is a general property of superconductors with broken 
particle- hole symmetry, also for temperatures far below To 

Defining tia = N' e A 2 /(AN e ), the total Hall force constant for a superconductor at u <C 
T< A becomes (see also the detailed discussion in Ref.i) 

1 + {uj T r y 

The first term describes the contribution arising from the quasiparticles bound to the core. 
Close to the superconducting-normal transition the scattering states have to be includedJl! 
and the term crosses over to the normal state Hall term; it therefore has the normal state sign. 
The second term is the hydrodynamic contribution. With uq = A 2 /ep and n&/n e ~ A 2 /e' F , 
the core term is dominant in the clean limit / > £(T), whereas the hydrodynamic contribu- 
tion determines the Hall angle in the dirty case. A sign change in the Hall effect occurs if 
the hydrodynamic term has a negative sign, i. e., N' e < 0. Within a free-electron based BCS 
theory we have N' e > and no sign change occurs. However, a simple modification of the 
electronic dispersion can drive N' e negative, resulting in sign changes of the Hall effect as de- 
scribed below (e. g., consider the dispersion = k 2 /2m + fc 4 /4m 2 eo in two dimensions: The 
corresponding density of states is A^ e (e) = m/ii(l + 2e/eo) and N' e = —(27i/meo)N 2 (e) < 
accordingly). With N' e < the Hall effect has the normal state sign in the clean limit and 
the opposite one in the dirty limit. Furthermore, the two contributions have a different 
temperature dependence through A(T), allowing for multiple sign changes. Our interpre- 
tation of the sign changes in HTSC (Bi- and Tl-based compounds) is as follows: at low 
temperatures the clean limit is realized with I > £(T) and the Hall effect has the normal 
state sign; with increasing temperature, r r and A decrease until the second term in Eq. @ 



dominates and a first sign change occurs when / ~ At even higher temperatures, close 

to Tq, the normal quasiparticles take over and a second sign change back to the normal 
state sign occurs. Note that the low-temperature sign change may be invisible if pinning 
is strong enough, which is probably the case for YBCO. This analysis provides a natural 
interpretation for the experimental findings as summarized by Hagen et a/.i. 

We now continue with an outline of the microscopic derivation of the vortex equation 
of motion, starting from a model Hamiltonian H that includes a short range attractive 
BCS interaction, as well as a long range repulsive Coulomb interaction (see Ref.@ for more 
details). We express the grand canonical partition function as an imaginary time path 
integral over the electronic fields tj) and the gauge field A a (a = r, x, y, z), 

Z = J V 2 ipVA a exp{-S} (3) 

with Euclidean action 

S = J d 3 r J dr(jp a [d T - ieA + f (V - ieA)}i/j a - 
-A-0 T ^;Vi^T + ieAorii + [E 2 + B 2 ]/8vr) . 

Here £(V) = — V 2 /2m — fi describes a single conduction band, and erii denotes the back- 
ground charge density of the ions. The idea is to construct an effective action for the vortex 
coordinate R only, by integrating out the electronic degrees of freedom. Our approach is 
inspired by the one of Simanek0. In addition to the analysis of Ref.0 we treat carefully 
the hydrodynamics of the problem and also avoid approximations for the matrix elements 
in the vortex core (see below). 

A Hubbard-Stratonovich transformation introduces the energy gap A as an order pa- 
rameter field and after performing the trace over the field ip (see RefsSlli for a survey of 
the technique used) we arrive at 

Z = J V 2 AVA a exp (Tr ln^ 1 - S ) , (4) 



Q- 1 



' d T - ieA + f (V - ieA) A ^ 
, A d T + ieA - f (V + ieA) ) 



So 



dx 




E 2 + B 2 

8^ 



+ ieriiA 



and J dx = Jq dr J d 3 r. The only remainder of the electrons is the Nambu-Gor'kov Green's 
function Q. The Euler-Lagrange equations obtained by varying Aq and A describe Thomas- 
Fermi and London screening respectively. They read0 



Both, the electronic density n e and current density j e are expressed through the electron 
Green's functions. For instance, n e = — Tr[a"3^] = \ J d£N e (£ + [/,+ ieA )[l — £/v / £ 2_ +~A 2 ~]. 
The electronic density is a function of the electro-chemical potential \x + ieA and in the 
presence of particle-hole asymmetry also of the energy gap Ailil. Due to charge neutrality 
n e (n, 0) = Hi. Expanding in Aq and A we find n e = rii + ieAoN e + ba + • • •, with wa = 
N' e A 2 /AN e . Deviations of n e from rii are screened on the Thomas-Fermi length \?f — 
(47re 2 A^ e )- 1 / 2 and yield a nonzero scalar potential Aq determined by the screened Poisson 
equation (— V 2 + \^p)A = Anien^. Magnetic fields and currents are screened on the scale 
of the London penetration depth = (47rn s e 2 /m)~ 1 / 2 . In the following we concentrate on 
strong Type II superconductors with \l 3> £• 

Varying A yields the BCS gap-equation, which has a constant as well as vortex solutions. 
Here we concentrate on the single vortex solution A(x) = A v (r — R(t)) with vortex coordi- 
nate R and A^ = | A„ | e l<fv . For A^ we adopt the mean- field solution from Ref .111. Using »a„ 
as a source, the screened Poisson equation defines also a single vortex solution A v0 for the 
scalar potential. In the limit of strong screening (Xtf <S A v0 = Ame\^ F n^ v . As a result, 
the electronic density in the vortex core does not differ from the density far away from the 
core. 

We neglect fluctuations around the mean field solutions A va {r — R(r)) and A„(r — R(r)), 
since longitudinal fluctuations of the phase and A are lifted to the plasma frequency, trans- 
verse fluctuations of A have a gap proportional to the superfluid density, and fluctuations 



V ■ E = 4iiie[n e (fi + ieA , A) — rii] > 



-d T E + V x B = 4yrj e . 



(5) 
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of |A| are at least at energy 2 Alia. Thus, the path-integral measure J V 2 AVA a reduces to 
JVR. 

Using a gauge transformation cAq — > eA$ — tp/2 = Qq and eA — > eA — V<£>/2 = 
Q, the energy gap A can be chosen real and manifestly gauge invariant quantities, such 
as the superfluid velocity Q/m, appear in £/ _1 . The dynamics of a vortex can now be 
studied by expanding TrlnC? -1 to second order in Q va and in the vortex displacement 
SA = — R(t) ■ VA„(r) around the static vortex solution. Furthermore, due to the singular 
gauge transformation a source term V x Vip v /2e = §o5(r ~ R) appears in the London 
equation that determines the magnetic field B„ around a vortex. 

We express the unperturbed Nambu-Gor'kov Green's function in the presence of one 
vortex in Bogoliubov-de Gennes eigenstates U\ with energy E\ as 



Ux(r)U{(r>) 



(6) 



In the relaxation time approximation that we use, uj^ = uj^ + sign(u; M )/2r r , where the uj^ are 
Fermionic Matsubara frequencies. 

The result of the expansion is an effective action S^iffR] = £>c + £>h for the coordinate 
R of one vortex, consisting of a hydrodynamic part 



S 



H 



+ 1 

dx( V Q _ — - + -HapQvaQvP — in>A v QvO j ■ 



8vr 



(7) 



and a core part 



S, 



c 



2(3- 



EE 



(R^W AV )(R_^W 



X'Xj 



XX' )J.,n \ luJ [J. 



Ex){iCo^ +n + Ex 



W 



XX' 



d 2 rUl 



VA,, 
VA„ 



Uy 



(8) 



and the u n denote Bosonic Matsubara frequencies. 

First we discuss the hydrodynamic contribution Sh- The kernel H a p is the polarization 
bubble and describes both longitudinal and transverse screeningil. The transverse part of 



the polarization term in Eq. (|7|) is 11^ = n s (T)(d~ ab — k a k b /k 2 ) + • • •, where a,b = x,y, z and 
the ■ • ■ indicate higher order terms in k and uj. The longitudinal part of the polarization 
term in Eq. ([?[) is ITl = N e + • ■ -. Using the equation of motion (|5|) for the fields and 
n e (n, A v ) — rii = ha„, the hydrodynamic part of the effective action can be rewritten as 

S H = -L [ dx (E 2 v + B 2 + X 2 L (-d T E v + V x B,) 2 + 

o7T J v 

X 2 TF (V ■ E v ) 2 + 16n 2 X 2 TF n 2 Av + Anin Av (p v ) . (9) 



The last term is the most important one for our discussion as it yields the hydrodynamic 
contribution to the Hall coefficient 7^ = 7rnA = ^N^A 2 /(AN e ). This result coincides with 
that of the TDGL-approach. The only other dynamic term in Eq. @ is the transverse part of 
the E 2 term that yields a small electromagnetic mas ME. All other terms are non-dynamic 
and contribute to the line energy of a vortex0. 

We now turn to the core contribution Sq in Eq. (||). Its origin is found in the transitions 
induced by the moving vortex between the CdGM states in the core labeled by A and A' and 
involving the matrix elements W\y. The energies are E\ = Xujq with A = half-integer. The 
sums over states A and A' in Eq. (§) may be evaluated using the constant level separation 
E\ — E\-i = cujii and properties of the matrix elements^. Explicitly we use the identity 



Uy = {Ey - E X )U{VUy , (10) 



VA,, 
VA,, 

together with the relations for the eigenstate wavefunctions V X U\ = (k F /2)[U\-i — ^a+i]> 

Ref.0 for a discussion of this point. Using the orthog- 
onality relations, we find the selection rule that only neighboring states are connected by 
the matrix elements. We also restrict ourselves to the temperature range ujq T where the 
sum over A's may be replaced by the integral / cIE\/ujq. The result for Sc can be written in 
the form 



S c = -[ dr j dr'\K+(r-r')R(r)-R(r') + 
2 Jo Jo L 

+iK c {t -t')z- (R(r) x R(t'))1 . (11) 



In Fourier components the kernels Kq determining the mass and damping (K + ) and Hall 
force (K~) are 

(12) 



_iuj n — ujQ zu n + uj _ 

where u n = u n + r r 7 1 sign(o; ri ). They are non-local in time, however, after analytic 
continuation to real frequencies they can be expanded in uj/ujq. The kernel Kq « 
— iu n (kpUQT 2 /2)/[l + (u;oT r ) 2 ], yields the core contribution 7 C as quoted in Eq. ([[]), if we 
put kp = 2im e in two dimensions. The kernel Kq rs (|o;„| + r r u; 2 )(A;f.a;oT r /2)/[l + (o;or r ) 2 ] 
is proportional to |a> n | for small frequencies, thus describing Ohmic dissipationHS Apart 
from the damping coefficient r/c it yields the core contribution to the vortex mas 

Mc = ^^mV (13) 



which is large in the superclean limit with ujqt t ^> 1. 

Thus, a complete description of intrinsic vortex properties can be obtained if both core 
and hydrodynamic contributions are included. The hydrodynamic part of the Hall-force 
was neglected in RefsJi3@0, whereas the core physics cannot be described by a hydrody- 
namic theory such as TDGlJ. In Ref.il the superconducting phase was coupled to the 
superfluid density in order to obtain a Galilei-invariant Magnus force Fm = k(\t — V) x z 
from hydrodynamics only. Our analysis shows that the phase of the superconducting order 
parameter couples to the square of the order parameter, with a small coefficient that de- 
pends on particle-hole asymmetry, i. e., details of the electronic band structure are relevant. 
A Galilei invariant Magnus force at T = and r r = oo in Fermionic superconductors is 
provided by the vortex core rather than the hydrodynamic flow around the vortex. 

In conclusion we have presented a microscopic derivation of the equation of motion of a 
vortex in a superconductor. Our results relate the observed sign change in the Hall effect 
in superconductors with I ~ £ to broken particle-hole symmetry in the electronic band 
structure. 
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